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Abstract 

We study the dynamics of brane worlds coupled to a scalar field and gravity, and find 
that self-tuning of the cosmological constant is generic in theories with at most two branes 
or a single brane with orbifold boundary conditions. We demonstrate that singularities 
are generic in the self-tuned solutions compatible with localized gravity on the brane: we 
show that localized gravity with an infinitely large extra dimension is only consistent with 
particular fine-tuned values of the brane tension. The number of allowed brane tension 
values is related to the number of negative stationary points of the scalar bulk potential 
and, in the case of an oscillatory potential, the brane tension for which gravity is localized 
without singularities is quantized. We also examine a resolution of the singularities, and 
find that fine-tuning is generically re-introduced at the singularities in order to retain a 
static solution. However, we speculate that the presence of additional fields may restore 
self-tuning. 



*J. Robert Oppenheimer fellow. 



1 Introduction 



There has been renewed interest in Kaluza-Klein theories over the past two years, mainly due 
to the reahzation that locahzation of matter ^ and locahzation of gravity Q may drastically 
change the commonly assumed properties of such models. These theories clearly open new 
approaches to the cosmological constant problem^, since using extra dimensions the no-go the- 
orem of Weinberg for adjustment mechanisms Q may be circumvented. A particularly simple 
scenario has been recently suggested in to at least improve on the cosmological constant 
problem (see also for an earlier mechanism involving extra dimensions), by proposing a 
mechanism for the cancellation of all order Standard Model (SM) loop contributions and the 
leading (tree- level) gravity contribution to the effective 4D cosmological constant, 

OjM^l + OjThrl + 0{TlM-^) + . . (1.1) 

tree-level SM loops higher order 

where M is the fundamental scale of gravity, the Planck scale in the bulk for instance. The 
mechanism is based on a single 3-brane (to which the SM fields are localized) embedded into 
a A+1 dimensional space-time. The essential new ingredient is a bulk scalar field (j)^ which is 
coupled to the brane tension T^^ (assumed to be the full loop-corrected SM vacuum energy 
density). The authors of [Q,^ then showed that one can find static solutions to the classical 
equations of motion for a vanishing bulk potential of the scalar field, but for arbitrary values 
of the brane tension Tft^. However, all the solutions found in [^,|5[ which localize gravity 
involve a naked space-time singularity, which has been interpreted as the boundary of the 
extra dimension (see also Q for a discussion on singularity in a brane- world context). Since 
the bulk is effectively compactified, we have to worry whether the size of the extra dimension 
is compatible with our phenomenological knowledge of gravity which has been tested up to 
millimeter distances. The proper distance from the brane to the singularity is found to be 

— 9 —1 —2/3 

He ~ Kg Tjjj, where M5 = is the 5D Planck scale, which is related to the 4D Planck 

scale, M4 = kJ^, by T^r = i^l/i^i- If the tension, associated to the vacuum energy of the SM 
fields, is of the order of the electroweak scale, Tf,r ~ TeV^, we naturally obtain]^ 

Ms ~ lO^GeV and yc ~ Imm (1.2) 

which is phenomenologically safe. Finally, while the SM contribution to the 4D cosmological 
constant would be of the order of 10~^^M4, the self-tuning mechanism cancels this term. 
However, the next possible term is of the order of 10~^^M4, which is still forty orders of 
magnitude too large. It is worth noticing that the self-tuning mechanism eliminates twenty 
orders of magnitude and thus should be considered on equal footing to the Randall-Sundrum 
solution to the hierarchy problem. 

The necessity of singularities in the bulk is here only dictated by the phenomenological 
requirement of localized gravity. However it complements the singularity theorem of Hawking 

*The cosmological constant problem concerns the brane cosmological constant that governs the expansion of 
our universe and it has to be distinguished to the vacuum energy density, or tension, on the brane: in particular 
it is possible, for a non-vanishing tension, to find solutions to the Einstein equations with a flat Minkowski 
metric on the brane in which case the brane cosmological constant is zero. 

^The fact that the electroweak scale on the brane is five orders of magnitude below the fundamental scale, 
Afs, is the gauge hierarchy problem in this model and a mechanism, such as low energy supersymmetry for 
instance, is needed to cancel the quadratic divergences which should bring Tbr near M^. 



1 



and Penrose under which generic initial conditions lead to singular solutions of Einstein's 
equations p. 

Several works have studied the self-tuning mechanism. In this paper, we examine 

the general properties of the solutions to the coupled 5D scalar-gravity system in detail. First, 
we investigate the solutions in the presence of a general bulk potential for the scalar field. The 
motivation to include a scalar bulk potential are twofold: (i) to give a mass to the scalar field 
and thus evade the cosmological problems associated to a massless scalar field coupled to gravity 
and/or those associated to an unstabilized extra-dimension; (ii) to overcome the singularity 
problem by considering a more general bulk geometry. We explain that the self-tuning behavior 
is expected to occur for a generic bulk potential. However, in the generic solutions (except for 
a vanishing bulk potential) self-tuning does have the more restricted interpretation that there 
is at least a finite region for the brane tension for which static solutions can be found (but not 
necessarily for all values of T^r). More precisely. Standard Model corrections are expected to 
occur at the weak scale, whereas self-tuning works up to the larger 5D Planck scale. 

After presenting our general results and methods on the solutions in the presence of a generic 
bulk potential we ask the following question: are there bulk potentials such that the self-tuning 
solutions avoid naked space-time singularities for a range of values of the brane tension in such a 
way, that gravity is localized to the brane (so as to reproduce 4D gravity for the observer on the 
brane). We show that for the 5D scalar-gravity system the only possibility for such solutions 
is when space-time asymptotes to five dimensional anti-de Sitter space (AdSs) far from the 
brane (therefore producing solutions of the sort considered in [|l^]). After a careful analysis we 
show that such solutions are always fine-tuned; that is, they occur only for particular isolated 
values of the brane tension Ti,r- We show that the number of allowed brane tensions is related 
to the number of negative stationary points of the bulk potential, which appears as a maximal 
index. In the case of an oscillatory bulk potential, we thus obtain a quantization of the brane 
tension as in a brane world construction from supergravity |12]. 

The paper is organized as follows: in Section 2 we present our method for searching for 
solutions of the coupled scalar-gravity system using the superpotential formalism of p, 13 T^. 
The advantage of this method is that it results in first order ordinary differential equations, 
and all the degrees of freedom are transparent without having to make a particular ansatz 
for the solution. In Section 3 we present a class of exactly solvable systems, given by an 
exponential bulk potential, which includes all the models presented in [^,|5|. We analyze these 
models in detail and find that all of these solutions necessarily involve a naked singularity or 
lead to an infinite Planck scale on the brane. We also find that, by relaxing an ansatz for 
solutions made in js), the exponential bulk potential is in fact self-tuning. In Section 4 we 
compare a perturbative and a numerical method to solve the equations for the most general 
bulk potentials, and present an example for both methods using the exactly solvable cases of 
Section 3. In Section 5 we present a no-go theorem for self-tuning branes that would localize 
gravity without singularities. In Section 6 we comment on the effect of resolution of the 
singularities on self-tuning. We conclude in Section 7. 



2 Self-tuning and scalar fields 

In this section we demonstrate that the 4D cosmological constant of branes coupled to scalar 
fields is generically self-tuned to zero, as in [§,|5|. We begin with the action for a brane coupled 
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to gravity and a real scalar field (such dilatonic domain walls have been extensively studied 
from a supersymmetric point of view [|l6|, a new ingredient here is the localization of SM 
fields)[|, 

S = J d^x^^\[R-%^^{dM<Pd''^ + V[<P])] -J d^-'x^l^J^MT, (2.1) 

where Gmn is the D-dimensional metric and g^i, is the pullback of the metric onto the flat 
domain wall at x^^^ = y = 0. For now we will not be concerned with the origin of the 
coupling f[(j)] to the brane, and allow it to be arbitrary; T is the non-canonically normalized 
brane tension and includes standard model vacuum contributions. We look for static solutions 
with the metric ansatzQ, 

ds^ = e-2^(j')/(^-i)dxl,_i + dy\ (2.2) 

The unconventional normalization of the action ( p.ip and the warp factor in ( p. 21 ) will be 
convenient in what follows. 

The equations of motion which follow from the action (^]^) with the ansatz (^]^) are [ pA| , [T5[ | , 

A"iy) = cl>'{yf + f[<Piy)]T6iy), (2.3) 
A'{yf = cP'{yf-V[^{y)], (2.4) 

ny)-A'm\y) = l'-^ + '-MT5iy), (2.5) 

where the primes denote derivatives with respect to y. If (/'(y) is monotonic in the bulk then 
the bulk equations of motion can be written in a first order form [11, 13, 1^] introducing an 
auxiliary field 

A'{y) = W[(^{y)]. (2.8) 

Because of the relation ( |2.6| ), W[(j)] will be called superpotential even though no supersymmetry 
is involved. We will use the first order formalism in order to construct solutions in the bulk 
on both sides of the brane: W±, (j)± and A± on the right (+) and left (— ) hand side. But 
the equations of motion (|2.6| )- (|2.8| ) must then be supplemented by boundary conditions at the 
brane. The boundary conditions due to the delta-function terms in ( p.3| ) and (2.5) can be 
written, 

_^ dW[cPiy)] _^df[cPiy)] 



d(t>{y) d<p{y) 
AA' = AW[<l>{y)]=Tf[cPiy)] 



{2S 

y=0 



*Our conventions correspond to a mostly positive Lorentzian signature ( — h . . • +) and the definition of the 
curvature in terms of the metric is such that a Euclidean sphere has positive curvature. Bulk indices will be 
denoted by capital Latin indices and brane indices by Greek indices. 

^We will not look for non-flat solution on the brane such as de Sitter or anti-de Sitter 4D configurations. 
When they exist simultaneously with flat Minkowskian solutions, it is a dynamical question to know which 
solution is preferred by stability. A nice feature of the case with vanishing bulk potential is that the flat 
solutions are the only 4D maximally symmetric solutions M]. 
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where AF indicates the jump of a discontinuous function F at y = 0, F{0^) — F(0^). In 
addition, (j){y) and A{y) must be continuous across the boundary. Hence, there are four 
boundary conditions at the brane. 

A count of free parameters in the solutions to the equations of motion immediately demon- 
strates that given f[4>{y)], the tension T is generically not fine tuned. If we do not impose 
orbifold boundary conditions in addition to those above, then there are naively six free parame- 



ters: one from the solution on each side of the brane of each of equations (|2.6D -(2.8). However, 
one overall constant shift in A{y) is not relevant because A{y) enters into the equations of 
motion only through its derivatives. That leaves five free parameters and four boundary con- 
ditions. There is generically a (finite) line of solutions for a region of scalar-brane couplings 

Once again, this is what we mean by self-tuning: Given an arbitrary scalar-brane coupling 
f[(j)] (possibly satisfying some constraints), there is a range of "brane tensions" T such that 
static solutions exist. In the generic case, as argued above, there is in fact a continuous set 
of static solutions for a given boundary condition specified by f[(p] and T. Furthermore, if 
/[0] = /'[(/>] then the range of T often extends to infinity. The reason is that if W[(p] ^ Vi'P] 
asymptotically when W is large, then W ~ W there, and T can be chosen arbitrarily large 
such that ifTJ'T) ~ {2W,2W'). 

Orbifold boundary conditions are more constraining. The additional constraints from the 
orbifold condition are A{y) = A{—y) and (p{y) = 4>{—y), which implies that VF+[(/)] = — 
where W± [(j)] are the solutions for W on the two sides of the brane. However, continuity of A{y) 
and (/>(y) is then guaranteed. Hence, in this case there are two free parameters (there would 
be three but a constant shift in A is not relevant) and two boundary conditions, and there is 
generically a solution for a region of couplings f[(f)] T. Thus, there is generically self-tuning in 
this case, as well. 

In the absence of orbifold type boundary conditions, if there are A'^ branes then there are 
3A^ -|- 3 — 1 = 3A^ -|- 2 free parameters and 4:N constraints, leaving 2 — N free parameters in the 
solution for a given set of boundary conditions. Hence, there can be up to two branes without 
fine-tuning. With orbifold boundary conditions, where image branes are included in N and 
there is assumed to be a brane at the orbifold fixed point (which contributes 1 to N), there 
are 3(iV - l)/2 3 - 1 = (3iV + l)/2 free parameters and 4(A^ - l)/2 + 2 = 2N constraints, 
leaving (1 — N)/2 free parameters in the solution. Hence, only if there is a single brane at the 
orbifold fixed point will self-tuning occur. If the space is compactified on a circle with orbifold 
boundary conditions, then a parameter count for the case of branes at the two orbifold fixed 
points demonstrates that a fine-tuning is necessary in this well Namely, there 

are two free parameters in the solution but four boundary conditions. In what follows we will 
concentrate on the case of a single brane coupled to a scalar field. 



3 Integrable bulk potentials 

In this section we study some special cases which were also partially discussed in Our 
discussion of the exact solutions with a vanishing bulk potential is similar, but we will not 
restrict ourselves to the ansatz A'[(f)\ oc cj)' made in |5|. In agreement with the parameter count 
in the previous section we will find that there is no fine-tuning in these theories, although some 
of the exact solutions exhibit non-generic behavior. 



4 



Let us first find the exactly solvable models. The challenge is finding a class of solutions 
to the nonlinear equation (|2.6D for a given V [(/)]. If (in a region wheref] y < 0) we write W[(/)] 
and W'[(l)] as|, 

w = iy3y[^(^Mc/>] + -L.), (3.1) 

w = iy^(M<A]-^), (3.2) 



then ( ^.6| ) is immediately satisfied for any prepotential 'w[(f)]. The consistency of ( |3.lD and (3.2) 
then translates into a differential equation for w[(j)]: 

Exact solutions for w[(f)] can be found when ( p.3| ) is separable, i.e. when is a 

constant. 

We distinguish three cases: 

• a vanishing bulk potential: V[(l)] = 0. 

• a negative bulk cosmological constant: V[(f)] = A < 0. 

• an exponential bulk potential: [(/)] = const. / 0. 

3.1 Vanishing bulk potential 

This case has been extensively studied by refs. Q|5| but it is a worthwhile exercise to repeat 
the discussion in terms of a superpotential. The equation for the superpotential can be solved 
without introducing a prepotential. Indeed, eq. (^) becomes simply 

W'[(^f = W[(^f, (3.4) 

with two branches of solutions, 

W[(l)]=ce"''. (3.5) 

where c is a constant of integration and e is a sign, both of them can take different values on 
the two sides of the brane (the constants of integration relative to the right (left) hand side of 
the brane will be denoted with a + (— ) subscript). With this form of the superpotential, the 
eqs. ( |2.7D -(^^ can be easily solved as 

•/•(y) = -eln(d- cy) (3.6) 
A{y) = -ln{d-cy) + e (3.7) 



The ca se o f a po sitive bulk potential can be studied in a very similar way up to some changes of sign in the 



equations (3.1)-(3.2) 



^From now, we will denote by a prime a derivative of V, f,WoTuj with respect to (p; or a derivative of ( 
or A with respect to y. 
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where d and e are some constants of integration that can also differ on the sides of the branes. 
Moreover, to make sense eq. ( ^.61 ) need: d^ > and d- > 0. Thus the continuity requires 

(j)(0) = (j)Q = -e+ In d+ = -e_ In d_ (3.8) 

A{0) = Ao = e+-lnd+ = e-- In d_ (3.9) 

while the jump equations are 

'-^-'-^=nMT (3.10) 

^-^^ = /MT (3.11) 

From the expression of the warp factor, we conclude that the Planck scale on the brane, 

-2 _ -2 

brane i.e. c+ > and c_ < 



1 ~ '"D I dy e 3)A/(D 1)^ jg f^j^ii^g ijff singularities are encountered on both sides of the 



• if e+e_ = 1: the consistency of the jump equations requires that /'[</'o] = e/i^Ao] and 
then it is possible to find solutions with or without singularities for any value of the brane 
tension but the singular solutions correspond to f[4)o\ T > 0. 

• if e_|_e_ = —1: the solutions with singularities exist only if /[(/)o]T > and —1 < 
f'[4>o]/f[4>o] < 1 but do not require any fine-tuning. 

If we choose f[(j)] = Ce'^'^ for the case e = e+ = e_ then it is clear that the boundary 
conditions can be satisfied for any T. There are several important comments to make about 
this case, which is quite non-generic. First of all, the fact that a specific form had to be 
chosen for f[(p] is a result of two non-generic features of the model: First, there is a symmetry 
(j) ^ (p + const.. As a result of this symmetry, one of the free parameters, namely (/>o, does 
not appear on the left hand side of the boundary conditions for the derivatives (f)' and A' . In 
addition, it turned out that the left hand sides of the two boundary conditions had the same 
form up to a sign. As a result, only f'[(l)]/f[(l)] is relevant, and given one solution {f[(p]T, f'[(p]T), 
there is an infinite set of solutions with the same f[(p] and arbitrary T. The fact that T is 
completely arbitrary is most likely unique to the case of vanishing bulk potential. But the 
fact that given /, /' ~ OiM^), self-tuning occurs for T to within ©(Ms) > 0{Mew) (where 
0{Mew) is the expected Standard Model contribution to the tension) is generic, and is what 
we mean by self-tuning. 

Even given the constraint on /[</>] from the shift symmetry in (j) in this case, as explained 
in IQ,^ the required exponential form of f[<j)\ might be natural from a stringy perspective where 
(j) is interpreted as the dilaton. 

Furthermore, as pointed out in j^, the solutions with singularities on either side of the brane 
must be chosen in order to have a finite gravitational coupling (assuming that the spacetime 
can be cutoff at the singularity in a consistent way). This feature will turn out to be generic 
except in theories which admit solitonic solutions for fine-tuned boundary conditions i.e. for 
special discrete values of the brane tension. 

3.2 Bulk cosmological constant 

Even though the case of a (negative) bulk cosmological constant can be studied without in- 
troducing a prepotential, we will present our method on this rather simple example before 
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proceeding to the somewhat more comphcated case of an exponential bulk potential in the 
next subsection. We will denote by V[4>] = A < the cosmological constant. The equations of 
motion are: 

^ = u, (3.12) 



^ 1 ,/^, ,n _ .,-2 

dA 



dy 



2v Atuil-LO-^), (3.13) 



— - ^V^uj{l + LO-^) . (3.14) 
dy 

The differential equation for the prepotential can be easily integrated, 

uj = ce't', (3.15) 

where c is a constant of integration. Plugging this expression for the prepotential into the 
remaining equations of motion, we obtain, 

^{y) = -Hct9{y)) , (3.16) 
A{y) = -ln\^{y)\+ln\l-^\y)\+a , (3.17) 

where a is a constant of integration and another constant of integration, y^, also appears in 
the expression of the function i?(y) defined, on the two different branches of solutions, by 

i?(y) = -tanh^(j/-ye) or ^y) = - coth ^ {y - y^) . (3.18) 

The nature of the solution depends on the sign of yc on each side of the brane. On the right 

(left) hand side, a positive (negative) value of yc will correspond to a solution involving a 
singularity at a finite proper distance from the brane, y = y^.- Near this singularity, the warp 
factor, exp{—2A/{D — 1)), goes to zero so the singularity appears as an horizon in the bulk. 
Conversely, if y^, the value of yc on the right hand side of the brane, is negative (respectively, 
y~ > 0), then the transverse dimension will be infinitely large; however, near infinity, the warp 
factor blows up and the Planck scale on the brane diverges, ruining any phenomenological 
relevance. 

The values of the constants of integration are constrained by the continuity and jump 
equations (with ■d± = tanh \/— Ai/^/2), 

(c_i?_) , (3.19) 
+ a+ = - \n\^-\ + ln\l - -dll + a- , (3.20) 

= f[MT, (3.21) 
^=f'[MT. (3.22) 

A solution to these equations can be found, provided that /'[0o] 7^ =t/[0o], for any value of 
the brane tension such that 



(t>0 




ln(c+^+) = -: 


Ao 




ln|i?+| +ln|l ■ 


1 

2 ^ 




1?+ i?_ 


1 

2 ^ 







7 



Moreover, singularities will exist on both sides of the brane iff 

fmT>0 and -1<1^<1, (3.24) 

fm\ 

just as in the case of a vanishing bulk potential. 

It is interesting to look at the Z2 symmetric solution for which there are the additional 
constraints: 

yt = —Uc ^ c+ = —c^ and a+ = a_. (3.25) 

The continuity conditions are automatically satisfied but the consistency of the jump equations 
requires a fine-tuning, 

fm'-f'm' = -^- (3.26) 

As already discussed in the case with a vanishing bulk potential, this fine-tuning is a conse- 
quence of the translational symmetry, (p ^ (p+ const., in the theory. As before, because of this 
shift symmetry we lose the appearance of a free parameter to adjust in the jump equations, 
which leads to a more restricted set of boundary conditions than for a generic bulk potential 
V[(f>]. The fine-tuning ( |3.26D is precisely the one appearing in the Randall-Sundrum model 
when the scalar coupling to the brane is a constant: 

^bk = ~ Qi^jj _ 2^ 4^feri (3.27) 



where A^fc and Thr are the canonically normalized quantities [12|, 



The solution constructed by Randall and Sundrum that localizes gravity with an infinitely 
large extra dimension corresponds to a limit of the singular solution where the singularities 
are pushed to infinity i.e. {)± = ±1. The jump equations then require 

/'[0o]=O and f[MT = 2V^ (3.29) 

i.e. the coupling between the brane and the scalar field vanishes and the canonically normalized 
brane tension, T^,., is fine-tuned to (3.27). In this limit, the expressions for the scalar field and 



the warp factor simply become 

(/> = (/>o and A = V^\y\. (3.30) 
3.3 Exponential bulk potential 

The last case that can be solved analytically with our method involves an exponential potential 
for the scalar field in the bulk. We will concentrate on negative potential while the case of 
positive potential requires some minimal changes. So the bulk potential will be parametrized 
by two real numbers a and b: 

V[^] = -a'^e^'"^ . (3.31) 
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On each side of the brane, the equations of motion are simply: 

duj , c<j^ + 1 

dy 



— = oo-bio— , (3.32) 

iae'"^a;(l-a;~2) ^ (3 33) 



rl A 

^ = iae^^a;(l + a;-2) . (3.34) 
dy 

The sign of a is not fixed and can be chosen independently on the two sides of the brane, as 
we will discuss. The first differential equation can be easily solved to express (j) in terms of oo: 

e-^* = e-'"^ |.;r^/(i+^) |l + 6 - (1 - h)u:''f/^'^'-^'^ , (3.35) 

where c is a constant of integration. This last result can be used to obtain a parametric 
representation of A and y as functions of w: 

, , . -.c r,- i^r^/^^^^^|i+6-(i-5)^t^^"^-^^ ,3 3,, 

yM = -.e jju: (1 + 5- (1-.)^^) ' ^'-''^ 

A{u) = A,- r dC. Jl^^P • (3-37) 

^0 w (1 + - (1 - b)uj^) 

On the two sides of the brane, the parameter a can differ by a sign while the initial bound of 
integration, ujq, and the constant of integration, c, can take any values compatible with the 
continuity conditions. A Z2 symmetric solution will correspond to two different choices of sign 
for a but the same values for ujq and c. 

Different kinds of solutions can be obtained depending on the value of b and of the range 
of integration for the variable a): 

• 6 < — 1 : in that case, dy/duj has an integrable singularity at a; = +00 but a non- 
integrable singularity at a; = while the singularities of dA/dto are both non- integrable. 
So we can find solution with or without bulk singularity at finite proper distance. 



— Solutions without singularity will be given by (^)-(|337|) where on the right (left) 



hand side of the brane, a has to be chosen positive (negative) and is a negative 
initial bound of integration. The parameter lo will range from u)q to 0~. It is easy 
to find the asymptotic behavior of this solution for large \y\, i.e. u near 0~: 

A ~ -ln|y| , (3.38) 

|j/|~oo 

from where it becomes evident that the singularity free solutions do not localize 
gravity because the Planck scale on the brane diverges. 

Solutions with singularity will still be given by eqs. ( 3.3(j )-( 3.37 ) with a positive 



(negative) parameter a on the right (left) hand side of the brane. And the range of 
integration goes from a positive initial value, ujq, to +00. In that case, y reaches a 
finite value yc while the warp factor goes like: 

A \n\y-y,\, (3.39) 

y^Vc 
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which indicates that the singularity is a horizon where the metric on the brane 
vanishes. The behavior of the warp factor near the singularity insures that the 
Planck scale on the brane, = J dy e~^^~^^^/^^~^\ is finite. 

— 1 < 6 < 1 : in that case the singularities oi dy / duj ai uj = ^oo , uj = ity^ (1 + 6)/(l — 6) 
and w = are integrable while those appearing in dA/dco are non-integrable. All the 
solution are singular with a horizon or a curvature singularity. 

1 < 6 : this case is quite similar to the first case because the singularity of dy/duj at 
u; = is integrable but the singularity at w = +oo is non-integrable while the two 
singularities of dA/dco are both non-integrable. So we can construct solutions with or 
without singularity. 

— Solutions without singularity will be given by (|]3|)-(|]3^) where on the right (left) 
hand side of the brane, a has to be chosen negative (positive) and is a positive 
initial bound of integration. The parameter lo will range from ujq to -|-oo. It is easy 
to find the asymptotic behavior of this solution for large \y\, i.e. u near -|-oo: 

A ~ -ln|y| , (3.40) 

|j/|~oo 

from where, once again, it becomes evident that the singularity free solution does 
not localize gravity. 

— Solutions with singularity will still be given by eqs. ( 3.36| )-( [337 ) with a negative 



(positive) parameter a on the right (left) hand side of the brane. And the range of 
integration goes from a negative initial value, loq^ to 0~. In that case, y reaches a 
finite value yc while the warp factor goes like: 

A \n\y-y,\, (3.41) 

y^Vc 

which indicates that the singularity is a horizon where the metric on the brane 
vanishes. 

• 6 = ±1 : in these two cases, we can construct singularity free solutions with a blowing 
up warp factor as well as singular solutions with a horizon at a finite proper distance. 

Figure |l| illustrates the different types of solutions. 

The main result of the analysis of these integrable bulk potentials is that we find two kinds 
of solutions: (i) solutions with a horizon in the bulk at a finite proper distance from the brane 
and with a finite lower dimensional Planck scale; (ii) solutions without singularity at a finite 
proper distance but associated to a bulk geometry that decouples gravity on the brane. Both 
kinds of solutions do not require any fine-tuning and can be constructed for a range of brane 
tension, T. However it seems impossible to find a singularity free solution that localizes gravity 
on the brane, unless the brane tension is fine-tuned as in the Randall-Sundrum model. This 
point surely deserves further scrutiny which the next sections will be devoted to. 
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Singulai' solutions Non-singular solutions 



Figure 1: Shapes of the warp factor for singular and non-singular solutions to the equations of 
motion with an exponential bulk potential. We have drawn Z2 symmetric solutions for different 
values of the tension on the brane: when the tension increases, the horizon becomes closer and 
closer to the brane for the singular solutions while the warp factor, e~^'^/^'^~^\ blows up faster 
and faster for the the singularity-free solutions. In both cases, the jump conditions require 
that the scalar coupling to the brane satisfies: \{df / dcj)) / f\ < 1 on the brane. 



4 Perturbative and numerical methods 

In this section we give two methods for finding approximate solutions to the coupled equations 
of motion in the bulk that satisfy the boundary conditions at the brane. First, we present a 
perturbative method, which we then show to break down for the interesting case of perturbing 
around a solution with localized gravity. Then we give a systematic numerical method for 
solving the equations. 



4.1 A perturbative method 

As emphasized in the previous section, the key to finding the self-tuned solutions is to utilize 
the fact that one can transform the superpotential without changing the bulk potential for the 
scalar field. Thus, one has an additional degree of freedom if one picks a different superpotential 
function on the two sides of the brane, and [(/>]. Let us now assume that we have 

found a static solution to the equations of motion (j)o{y) and AQ{y), and assume that we have 
chosen one of the integration constants such that Aq{Q) = 0. This is a solution for a particular 
value T of the brane tension. In order to find the solution obtained by perturbing the brane 
tension as T — > T + 5T, we first need to find how one can change the superpotential W around 
VFo[(/>] so as to leave the bulk potential unchanged: 

v[<p\ = wl,[<pf-wM\ (4.1) 



which should be invariant under W ^ W + 6W. Linearizing (4T) around Wq we obtain a 
differential equation for SW of the form 



6W' W( 
^0 



5W ~ W'' ^^-^^ 
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which is solved by 



5W[ct>] = Cexp I (4.3) 



where the arbitrary constant C yields the extra degree of freedom needed to find a solution 
for any value of T. The superpotential variation can be expressed in terms of the unperturbed 
background solution using the equations of motion d<j)o/WQ[(j)o] = dy and VFoifAo] = dA^/dy, 
we obtain 

5WWy)] = C exp / dy'^ = C e^^'^^l (4.4) 
J dy 

From ( [4. 3D we can see that the change in the derivative of the superpotential is given by 

6W'[<t>{y)] = C^e''^^y\ (4.5) 

my) 

In order to satisfy the jump equations ( |2.9| ) for the perturbed tension T + 6T, we need to 
choose C differently on the two sides of the brane. The values of C± are then given by 

where AAQ/(f)Q denotes the jump in Aq/(Pq at y = 0. Once C± is determined from (|4.6| ), we 
can simply integrate the equations 

S(l)'{y) = 6W'[My)], 6A'{y) = 6W[My)] (4.7) 

to obtain the perturbed solutions (j)o{y) + 5</)(y) and Ao{y) + 6A{y). This method always 
results in a perturbed solution. However, in the most interesting case, when the unperturbed 
solution asymptotes to AdS space thereby localizing gravity to the brane (that is for A{y) ~ 
{D—l)\y\/ Ry^iiS for large values of y), one can easily see that the perturbative method presented 



here always breaks down. This can be seen by inspecting (4.4), which shows that in this case 
(5PF[(/)(y)] oc e^^""^^!^!/^^'*^ . Therefore the perturbed values of 6(f) and 6 A grow exponentially, 
and thus the linearized approximation breaks down. We will examine the case of localized 
gravity in detail in Section ^. But first we give a numerical method that can be used for any 
choice of the bulk potential. 



4.2 A numerical method 

Next we present a method for solving the system (p.6|)-( [2^ ) numerically for any brane tension 
and potential in the bulk. Thus the input functions are the bulk potential V[(p], the brane 
tension T, the coupling of the scalar to the brane determined by the function /[</>], and in 
addition we can pick the value of the scalar field at the brane 4>{0) = (f)Q arbitrarily. In order 
to find a numerical solution to these equations, one has to first make sure that the boundary 
conditions that one imposes do satisfy the jump equations (l2.9|) . Our strategy is the following: 
we first determine the superpotential functions W4.[0] and to the left and the right 
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of the brane numerically such that the boundary conditions arising from the coupling to the 
brane are satisfied. This can be done by noting that once (po is fixed, the jump equations are 
just given by 

W+-W^= foT, W'_^-W'_ = fl^T, (4.8) 

where W± refers to the values of the superpotential functions to the right and left of the 
brane at (po^ fo = f[4'o]j etc. In addition, the superpotential functions must be such that they 
reproduce the correct value of the bulk potential at the brane: 

W\.^ -Wl = Vo, W'J -W^= Vo, (4.9) 



where Vq = V[(f)o]. Eqs. (|4.8| ) and (|4.9| ) together are enough to determine the values of both 
W± and Wj_ at the branes. They are given by the expressions: 



Due to the quadratic nature of equations ( [4.§| ) and ( [4.S| ) there is a second solution, where the 
signs in front of the square roots in ( [4.10 ) are both simultaneously flipped. Once the value of 
W± and VF4 are fixed, one can numerically integrate the equation^ 



W^4[(^] = sgn(M^4)yFp+W±p2 (4.11) 

to obtain the superpotential functions to the left and the right of the brane that satisfy all 
boundary conditions. Once W±[(/)] are numerically known, we can simply integrate the equa- 
tions 

cP'iy) = W'[cP{y)], A'iy) = W[cP{y)] (4.12) 

to the left and the right of the brane to obtain the numerical solutions for 4'{y) and A{y). 
We will show an example for this below for the case when the exact solution is known, and 
compare the two results. 



4.3 An example for the perturbative method 

In this Section we test how well the perturbative method described in 4.1 works. We will 
compare the analytic solution of the model with a vanishing bulk potential to the perturbed 
solution around a different analytic solution. The main conclusions are as expected: the 
perturbative method works well far from the singularities. However it gets worse as we approach 
the singularity itself, and does not capture the essential feature of the self-tuning solutions: 
whereas in the self-tuning mechanism, for a fixed value of the scalar field on the brane, the 
place of the singularity adjusts itself with respect to the value of the brane tension, here the 

*We will see in Section 5 that, unless the brane tension is fine-tuned, the superpotential, W[<^], will be a 
monotonic function of (ji and thus W± [0] will keep the sign of W± . 
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singularity of the perturbed solution remains at the same place where the singularity of the 
unperturbed solution was. Therefore, we conclude that this method is not very efficient in 
capturing the basic properties of the self-tuning solutions. 

The example we consider is the vanishing bulk potential discussed in 3.1, with the choice 



e+ = l,e 
choose T 
given by (I 



T, and we choose f[(j)] = e' 



0/2 



1, and for the unperturbed solution we 



1, while we pick 6T 
-Q, with 



0.1 for the perturbed solution. The analytic solution is 



d± 



4 , c- 



(4.13) 



where A has been normalized to zero on the brane. The perturbed solution from (4.7) is given 
by 



T 



(4.14) 



where 5(j)± has been normalized to zero on the brane. The perturbed solution obtained for 
T = 1, 5T = 0.1 compared to the exact solution for T = 1.1 can be seen in Fig. ^. As mentioned 
above, the perturbative solution nicely follows the exact solution away from the singularity, but 
deviates from it close to the singularity, in particular the place of the singularity is incorrectly 
predicted to coincide with the singularity of the unperturbed solution. 




-2.5 



Figure 2: The exact solution versus the perturbed solution for the case of a vanishing bulk 
potential. The curve that blows up at smaller values of y corresponds to the exact solution. 
The initial unperturbed solution is given by the dashed curve. The singularity of the perturbed 
solution appears at the same distance from the brane as the singularity of the initial solution; 
the shift of the singularity with a variation of the tension is missing. 



4.4 An example for the numerical method 

We have shown above, that the perturbative method in general does not do a good job in finding 
the solutions, since it becomes unreliable close to the singularities. However, the numerical 
solution should not have these problems. Indeed, we analyze the same example as above (the 
case with vanishing bulk potential) using the numerical method, and find that the exact and 
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numerical curves are virtually indistinguishable. Therefore, we suggest that in order to analyze 
potentials for which no exact solutions can be found, one should use the numerical method 
rather than the method based on perturbations. 

We are looking for a numerical solution to the case analyzed perturbatively above, that is 
vanishing bulk potential, = e*^/^, <pQ = 1, T = 1 and e+ = l,e_. = —1. From Eqs. (3.4) 
and (B.SI) we find the starting values of the superpotential to the left and the right of the brane: 



4 ' 



4 



Numerically integrating the equation 



±W[(p] 



(4.15) 



(4.16) 



with the boundary conditions VF±[(/>o] = W± one obtains the numerical values for W±{(j)). 
Finally, the values for (j){y) can be obtained by numerically inverting the integral 



y 



00 



W'{(l)) 



(4.17) 



to the left and right of the brane. The numerical solution obtained this way overlayed on 
the exact solution of Section 2 can be seen in Fig. ^. One can see that the two curves are 
virtually indistinguishable, suggesting that the numerical method works very well around the 
singularities, and should be the preferred method of looking for solution in the absence of exact 
solutions. 



^(y) 

lo^ 
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Figure 3: The numerical solution overlayed on the exact solution for the case of a vanishing 
bulk potential. The fact that the two curves are indistinguishable shows that the numerical 
method works extremely well even close to the singularities. 



5 Localized gravity without singularities 
5.1 A no-go theorem 

We would like to reexamine in this section the count of free parameters versus fine-tuning 
parameters needed to preserve a static Poincare invariance on the brane, with the restriction 
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that the solution corresponds to an infinitely large extra-dimension (without singularities) in 
the bulk and localizes gravity on the brane. 

Consider the general D-dimensional background preserving Poincare^^i 

The graviton zero-mode is localized on the bran^ precisely when the effective Planck scale is 
finite on the brane |15|. In terms of the warp factor, this condition is: 

dy e-(^-3M(2^)/(^-i) < oo . (5.2) 



.2 ^2 



It is convenient to introduce a transverse coordinate z for which the bulk metric is conformally 
flat: 

ds"^ = Q'^{z) {dx%_i + dz^) , (5.3) 

where the conformal factor, is related to the warp factor, e~'^^^y^^^^^^^~^\ by the two 
identities: 

n{z) = e-^fe)/(^-i) and n\z)dz^ = dy^ . (5.4) 

Then the condition ( |5.2| ) is equivalent to having a massless normalizable bound state, which is 
interpreted as the graviton on the brane: 

V'o oc ^i^j^ J |^^|2 ^ ^ ^ ^5 5^ 

Let us assume that the behavior of ipo at infinity is a power law: 

V'o oc z"" . (5.6) 



The localization of gravity ( |5.5D then requires: a > 1/2. 

In our study, the value of the parameter a is constrained by the fact that the background 
is created by a scalar field coupled to gravity. From the equations of motion, we easily deduce 
that A has to satisfy: d'^A/dy'^ > 0, which translates, in the z coordinate, in a lower bound on 
the value of a: 

D-2 , , 

a > ^— . (5.7) 

Furthermore it is worth noticing that an upper bound on a comes by the requirement of 
a geometry without singularity at a finite proper distance^. Indeed the proper distance from 
the brane to infinity is given: loo = J dzQ that diverges iff: 

a < . (5.8) 



*We do not consider the recently proposed possibility that gravity might be quasi-localized to the brane 
since in those models the A" > condition is not satisfied, therefore it is not possible to generate those 
backgrounds from a single scalar field. 

^For a power law conformal factor, the curvature always vanishes at infinity. However quadratic invariants 
such as -RaiiM2M3M4-R*^^*^^*^^*^'' will be singular at infinity as soon as q > (_D — 2)/2. 
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So the only background for the scalar field coupled to gravity that localizes gravity without 
singularity (with an infinitely large extra dimension) is asymptotic, at infinity, to the horizon 
of an anti-de Sitter space, as in the RS model, and corresponds to a = (D — 2)/2. In that case, 
the warp factor is exponentially decreasing with the proper distance to the brane: 



A ^ {D-l)\y\/RAds ■ 



(5.9) 



The aim of this section is to show that such a background necessarily requires a fine-tuning 
between the brane and the bulk. This is not to say that there is no self-tuning in these models, 
only that the nonsingular solutions require fine-tuning. 

The previous asymptotic behavior has a nice interpretation in terms of the superpotential, 
[(/)]. According to the equations of motion, the fact that A is asymptotically linear means 
that (j) becomes constant and we will denote by (/>~ and (p'^ the asymptotic values of at 
y = —CO and y = +00 respectively. 

The equation 



d(t) _ dW 
dy d(j) ' 



(5.10) 



is similar to an RGE with W playing the role of the /3-function. In order for cp to approach a 
constant (fixed point) at infinity, the /3-function dW/d(j) must have zeroes. In other words, in 
order to extend the range of the transverse coordinate from y = — cxd to y = +00, the values 
<f)~ and <j)'^ at infinity must be some roots of dW/d(j): 



dW 



and 



dW 







(5.11) 



Furthermore, at infinity, A has to be linearly increasing in \y\; otherwise the conformal infinity 
of AdS would be reached without localized gravity |12|. Given that 



we conclude that cp^ and (p^ must satisfy, 



W[(p-] < and W[(pt] > 0. 



(5.12) 



(5.13) 



Finally, (p^ and (p^ must be dynamically reached at y = —00 and y = +00, which according 



to ( ^.10 ) is possible iff: 



d^W 



> and 



d^W 



<0, 



(5.14) 



or at least a similar condition for the first non-vanishing higher order derivatives at (p^ and 
Pictorially, the previous conditions are summarized in figure ^. 

Finally, the last equation of motion that relates the superpotential, W, to the scalar po- 
tential in the bulk, V, partially fixes the possible values of (p~ and <p'^. Indeed this differential 
equation evaluated at infinity gives: 



-w\ 



< and V[(pt] = -W[(ptf < 0. 



(5.15) 
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w = o 



w = o 



Figure 4: Asymptotic behaviors of W[(f)] leading to a singularity free bulk geometry localizing 
gravity on the brane. The absence of singularities is equivalent to the conditions: W ^ = 
W ^ = 0; the dynamics of the equations of motion require W ^ > and W ^ < 0, while the 
localization of gravity requires W~ < and > 0. 



while a differentiation with respect to (p gives: 
dV dW (d^W 



W 



thus 



dV 



and 



dV 



0. 



'\'t>t 



(5.16) 



More information on W can be obtained by considering higher order derivatives of the 
differential equation between W and V . Indeed it is easy to prove by induction the following 
relation: 



y W = ^ 2 ^\ W^^^ (^(n-k+2) _ (n-fe) 



(5.17) 



where y^"^ denotes the n*'^ order derivative of V and similarly for VF^"); in addition, we will 
denote by W^'^\ the values of VF^"^) a,i (j) = <; 



At the second order, by evaluating ( 5.17] ) at 



we obtain a quadratic equation for T^(^) 



The superpotential will be real-valued provided that: 

v^'-^t > K^/2 , 

and then there are four different branches at each asymptotic point: 



(5.18) 

(5.19) 

(5.20) 
(5.21) 



with ei = ±62 = ±1. 
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However the compatibility between the gravity locahzation requirement ( ^.130 and the 
dynamics of the differential equation ( ^.141) can be fulfilled iff 

> , (5.22) 

and only one out of the four branches is retained: 



l^+ = ^/^ Wi^ = ^ - . (5.24) 

At higher order, the relation ( |5.17| ) becomes linear in W^'^^^ and allows to compute VF^")^ 
recursively in terms of lower derivatives: 

n-l 

V("';^-^2(^:i)VK('=)^(H^("-'=+2)±-VK("-'=)^)+2{n~l)iy(2)±W'("-2)± 

Of course, for this expression to make sense for any integer n > 2, it is important that no 
solution to the equation nW^'^\ — = can be found. However, the roots of this equation 
would satisfy: 

which is incompatible with the physical requirements ( ^.15| ) and (5.22) of localized gravity. 

At this stage, it is worth noticing that the uniqueness of a superpotential W reaching a 
given asymptotic point follows from our requirements of a localized gravity without singularity. 
Otherwise, a scalar field V in the bulk could be constructed such that the equation nW^'^\ — 
Wf = has some solution in which case W^"'\ would not be determined, leading to a 
continuum of solutions. 

Finally, the whole expression of W reaching the asymptotic points cf)^ can be uniquely 
reconstructed from its derivatives through a Taylor expansion: 



oo ^ 

<0: 1^_[0] = ^ T^W;(</>-<^-)-, (5.27) 

n=0 
oo ^ 

>0: VF+[^] = ^ H^W+(</>-<^+)" . (5.28) 



n=0 

From the expression of the superpotential, we can solve the equations of motion on the two 
sides of the brane. And the jump conditions are given by: 

The first equation will fix the value, 0O) of the scalar field on the brane while the second one 
fine-tunes the value of the brane tension. Generically, we will obtain only discrete values of (pQ. 
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Indeed, if there exist a continuum interval of solutions for (pQ, then the jump equation becomes 
a differential equation that can be integrated on this interval and we obtain that f[(j)] has to 
be proportional to — VT- [<?!)], but in that case there is only one possible value for the 

brane tension^In all the other cases, given two asymptotic points, (j)^, we will obtain discrete 
solutions for (pQ and T. Moreover, using different critical asymptotic points, we usually find 
different values of the brane tension and the number of values of T that allows an infinitely 
large extra dimension with localized gravity is related to the number of critical points of the 
bulk potential V satisfying ( |5l5| ), ( pB] ) and ( [5:221 ): 



HT < i3nc - 2)ns, (5.30) 

where nj- stands for the number of values of T such that gravity is localized without a singu- 
larity, while nc is the number of critical points of V as defined above. The multiplicity factor 
comes from the fact that the scalar field can asymptote either the same critical point or two 
adjacent ones on the two sides of the brane. The upper bound has been semi-quantitatively 
corrected by taking into account the average number, ns, of solutions to the jump equation 
for (f)o. It may also happen that some values of T are degenerate. We will explicitly describe 
an example in the next subsection. 

It is worth noticing that in the case of an oscillatory bulk potential, we will obtain an 
infinite number of discrete values for the brane tension that is as quantized. 

To complete the proof of the no-go theorem presented above, we need to show that there 
is no loophole in the above argument due to the fact that we have used the superpotential 
formalism. The subtlety that one might worry about is that in the proof above we have 



implicitly assumed that (p is monotonic, by writing the second order equations ( |2.3D -(2.4) in 
terms of first order equations involving W. In particular, if 4> is not monotonic (that is if 0' = 
at a finite value of y) one does not have a globally defined superpotential function W{<p), but 
instead one must define separate superpotential functions Wi for the regions between yi and 
yi+i, where (p'iyi) = i?!>'(yi+i) = 0. For these superpotentials that are not globally defined it 
is then possible to have VF'((/)*) = without satisfying F'(0=k) = 0. However, it is impossible 
to continue the solution beyond This by itself however may not be a problem, as long as 
at (pif one is smoothly switching over to another branch of VF. Of course this switch-over can 
only happen at a point where <j)' = 0, since otherwise (j) is monotonic, and one can solve the 
equations in terms of the superpotential, which is well-defined around (p^. Next we show that 
such possibilities do not get around the no-go theorem presented above. The reason is that in 
order to have localized gravity with an infinitely large extra dimension, we need (f)' 0, cj)" — 
for \y\ oo, and therefore we find from the second order equation ( |2.5| ) that V' 0. Thus 
the critical points at infinity must belong to an "ordinary branch" described above, where W 
can be continued at both sides of (j)c- However, once we are on an "ordinary branch" which 
can be globally defined, all the critical points will actually happen at V' = 0. Since the only 
possibility for switching over to another branch is at W = 0, and at those points V = 0, 
one can never switch off the ordinary branch, and therefore one can not circumvent the no-go 
theorem by gluing non-monotonic (j)^s together. 

* Whereas the non-canonically normalized brane tension, T, is fine-tuned, the physical brane tension, 
Tbr tx f[(j}o]T, is not fine-tuned since the value of (j)o can vary continuously. Whether this is a solution to 
the cosmological constant problem or not, beside the fine-tuning requires on f[(j)], depends whether SM loops 
will modify T or Ttr- This issue deserves further analysis in future work. 
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5.2 A numerical example 

In this section we demonstrate the ideas of the previous sections in a numerical example. The 
bulk potential in this example (Fig. ^) is 



V{cP) = -(P^ + 110^ - + 1, 



(5.31) 



which is generated by the superpotential W[ 




Figure 5: The scalar bulk potential (5.31). The physically interesting region, near the two 
negative stationary points, is emphasized. 



As in some of our previous illustrative examples the potential ( |5.31| ) is unbounded from 
below, and so the theory might be unstable to quantum fluctuations. However, we will only 
be concerned with static solutions and for our purposes any instabilities will not be relevant. 

As emphasized in Figure ^, the bulk potential ( 5.31| ) has two negative stationary points 
at (/> = ibl/\/3. According to our no-go theorem, there are isolated superpotentials solving 
the equations of motion with critical points at ±1/^/3. Here these superpotentials are very 
simple because from ( 5.25| ) only a finite number of derivatives are non- vanishing thus the 
superpotentials are polynomial and take the form 



W[cP] 



+ 



reaches 4>c = =bl/\/3 on the negative (positive) branch, (5.32) 
reaches (/>c = ±l/\/3 on the positive (negative) branch. (5.33) 



Depending on which critical point we want to asymptote at infinity, we can construct four 
types of solutions that localizes gravity 



for 



1 / \/3, 0c =1/ V^) ■ the solution is 
Hy) = 75tanh(^/3|y-y±|) 



A{y) = ^tanh\V3{ 



y - yt)) + ilncosh(\/3(7/ 



yc )) + «± 



(5.34) 
(5.35) 



where and a± are four constants of integration to be determined by the continuity 
and jump conditions. The continuity conditions imply that = —y~ and a+ = a_. 
The jump equations will depend on the precise form of f [(/)]■ Except in the degenerate 
case where f[(f)] oc VF+[0] — that has been discussed in a footnote on page 2C, we 
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will generically obtain a finite number of discrete values for (pQ. For instance in the case 
of an exponential coupling, /[(/>] = ae*"^, the values of 0o will satisfy 



^o(</>g - 1) 



(5.36) 



which admits three solutions whatever the value of b is. And thus there exist three values 
for the brane tension that will lead to a localized gravity. 

for {(f)'^ = l/\/^,(f>~ = — l/-v/3) or {(f)'^ = —\/^/^^cj)~ = l/\/3) : in those cases, there is 
no discontinuity in the superpotential and the brane tension has to vanish. 

for (0+ = -1/^/3, = -l/\/3) : this case is analogous to the first one and, for an 
exponential coupling, three values for the brane tension are possible and they are just 
the opposite of the ones obtained in the first case. 



w[(M 




Figure 6: Solutions for ^(y) and A{y) in the bulk potential ( 5.31| ). There is a unique, 

fine-tuned, regular solution which approximates the Randall-Sundrum solution for large \y\. 
All other solutions are singular. 



Besides the previous fine-tuned solutions with an infinitely large extra dimension and local- 
ized gravity, all other solutions will either decouple gravity on the brane or involve a horizon at 
a finite distance. We would like now to examine numerically the self-tuning of these singular 
solutions. To do that we scan for solutions to (2.6) by fixing the value of W[4>\ at = 0. The 
solutions are plotted in Fig. Note that in agreement with the no-go theorem of the previous 
section there is a unique solution in this branch of solutions with stationary points at the 
positions of the minima of 1^ [</>]. As a result of the uniqueness of the solitonic solution, there 
are no solutions (in this branch) with < |1^[0]| < W^[</'*] — .385, where 0=,, is the value of the 
field at the position of the local minimum of V[(j)\, in this case 0=,, ~ —.577. Also in agreement 
with the results of previous sections, the solitonic solution is the unique regular solution for 
4>{y) and A{y), which in this case approximates the Randall-Sundrum solution for large \y\. 

In order to determine the range of boundary conditions at the brane which could be satisfied, 
we note that except for the solitonic solution, all other solutions span an infinite range in (/>. As 
noted in Section 4.2 the boundary conditions can be expressed in terms of W[(j)\ and M^'[(/>] at the 
boundary. Hence, if there is a space- filling range of solutions (VF [(/>], T^' [0]) then no fine-tuning 
is necessary in order to have a static solution with arbitrary boundary conditions in that range. 
In other words, a given (21^[(?!)o], 21^'[i?!>o]) can be equated with a certain (orbifold) boundary 
condition (/[i;^o] /'[c/'o] T). Then, given one such (/[c/'o] /'[(/>o] T) there is a continuous set 
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of solutions around that point. That means that for a given (/ [0o] ) /' ['?^o] ) i^i & certain range 
of /'// there is a range of tensions T for which there is a solution; hence the theory is self- 
tuning. Fig. 1^ is a parametric plot of W[(j)] versus for several solutions of It is 
intended to illustrate the fact that there is a space- filling region for W > .385. Given any 
boundary condition parametrized by f[(j)] and T, T can be rescaled by an amount given by 
the intersection of the set of solutions {2W, 2W') with a line through the origin and (/T, f'T). 
As is generically expected, one can see from Figure |^ that if /, /' are 0(1) = 0{M^), then 
T can be rescaled by 0{M^) without eliminating a solution. Hence we have demonstrated 
the self-tuning of this model. However, the caveat is that because of the isolated solitonic 
solution, there is a fine-tuned region near W, W ~ 0{1). If / is 0(1) at some matching point 
00, and T is 0{Mew/M^) ^ 1 then a fine tuning is reintroduced because of the uniqueness 
of the solitonic solution. There is still a large region of parameter space where the theory is 
self-tuning, but whether that region is natural or not requires exploration. This phenomenon 
is a result of the existence of solitonic solutions, and is nongeneric. Furthermore, if we do not 
require orbifold boundary conditions, then once again self-tuning is natural. Note also that 
asymptotically the parametric plots of W vs. W include the line W = W . This is generic in 
any region where the solutions satisfy VK[<^] ^ ^ [</•]• This also implies that for any solution 
/ ~ /', there is a very large range for T for which there are solutions. This behavior for large 
values of W is common, and extends the range of T over which self-tuning occurs possibly to 
iboo if f[(f)] = f'[(t>], which may be natural from a stringy perspective iQ,^. 

6 Resolution of singularities and fine-tuning 

In this section we reconsider the resolution of singularities^ as proposed in 
have seen, the case of zero bulk potential, which is the case studied in |l^, is quite non- 
generic. There is a shift symmetry in the scalar field which makes the boundary conditions 
more constraining, and two of the boundary conditions turn out to have the same form. Hence, 

* As they stand in our solutions, Einstein's equations are not satisfed at the singularity but require a singular 
stress-energy tensor located at the singularity that may correspond to a brane: the introduction of this brane, 
for instance, is what we mean by resolution of the singularity. 



IC, 17]. As we 
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we study the generic case, and propose a new resolution of the singularity which may restore 
self-tuning. 

The idea of |1^, 17 1 is to add a brane at each of the singularities such that the equations 
of motion, or boundary conditions, are satisfied there, as well. As pointed out in |10|, 17|^ 



the singularities contribute to the effective 4D energy density when integrated over the ex- 
tra dimension, and the contribution of the singularities is essential for vanishing of the 4D 
cosmological constant. 

However, addition of branes at the singularities adds new boundary conditions. In order 
to answer the question of whether or not the theory is self-tuning we must better understand 
the continuation past the singularities. In it is proposed that the spacetime is either 
periodically continued or cut-off at the singularity. In the case that the spacetime is cut-off 
at the singularities on each side, there are generically two new boundary conditions at each 
singularity, but no additional free parameters. Hence, without orbifold boundary conditions 
there are 3 + 3 — 1 = 5 free parameters and 4 + 2 + 2 = 8 boundary conditions, and the 
system is overconstrained. If one imposes orbifold boundary conditions, then there are only 
half as many additional boundary conditions (2 + 2 = 4 boundary conditions in all), but only 
3 — 1 = 2 free parameters. In either case the system is overconstrained. Hence, a fine-tuning is 
required. Although in the absence of a bulk potential the situation is modified because of the 
non-generic features mentioned above, in that case a fine-tuning at each of the singularities is 
required, as well. 

^^^"^ identified 





horizon horizon 



Figure 8: Penrose diagram for causally disconnected regions glued together at the resolved 
singularity. 

Although the details of any continuation of the spacetime beyond the singularity will depend 
on quantum gravitational dynamics, we propose another scenario which does slightly better 
than the cut-off scenario, although fine-tuning will still be required. The singularities are 
generically horizons (c./. Sections |3.3| ) because the warp factor e~^^(^)/(^~^) vanishes there. 
Hence, light does not cross the horizon and we can imagine a scenario in which there are 
bulk fields living beyond the singularity, out of causal contact with our universe. The Penrose 
diagram for this scenario is illustrated in Figure |8|. However, the bulk fields across the horizon 

^See also the first reference in for a discussion on vanishing 4D effective vacuum energy. 
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would provide an additional three free parameters which one might hope would help avoid 
fine-tuning. More precisely, following the counting in Section |2|, if there are orbifold boundary 
conditions, then there are 3 + 3 — 1 = 5 free parameters from the bulk fields on both sides 
of the singularity, but 2 + 4 = 6 boundary conditions. Hence, although the system is less 
constrained than the case in which the spacetime is cut-off at the singularity, one fine-tuning 
is still required. 

One might suspect that additional fields would add additional degrees of freedom which 
could help in self-tuning. For example, the addition of a scalar field with second derivatives 
in its equation of motion would contribute two additional free parameters on each side of a 
brane, but only two boundary conditions (continuity and change in the derivative at the brane). 
Hence in a system of branes there are net, at least naively, two additional free parameters from 
the extra scalar field, which could be used to restore self-tuning at the singularities or perhaps 
even produce nonsingular solutions. However, a more detailed analysis is required in this case. 



7 Conclusions 

We have studied brane worlds coupled to a scalar field and have found that self-tuning is a 
generic feature of these models. In these models the dynamics of the scalar field provides 
additional degrees of freedom, which generically alleviates the need for fine-tuning of static 
solutions. We have reexamined the exactly solvable models, two of which were studied previ- 
ously [@,|5|, and have found that those case are more constrained than the generic case because 
of a shift symmetry in the scalar field in these models. Still, these theories are self-tuning, 
including the case of an exponential potential. Whereas in a fine-tuning was necessary in 
this case due to a particular ansatz for the scalar and graviton fields, we showed that the more 
general solution is not fine-tuned, in agreement with a counting of free parameters in these 
models. We demonstrated that singularities in the self-tuned solutions are generic if gravity 
is to be localized, and we presented a no-go theorem to this effect. We provided perturbative 
and numerical techniques in order to calculate the self-tuned solutions, and we illustrated the 
major points of the paper via several numerical examples. Finally, we pointed out that the 



fine-tuning that is required in order to resolve the singularities in the spirit of |10|, 17| for the 
case of zero bulk potential is generic. 
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